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i. Basis of the report 

1 . This report has been drawn on the basis of (substitute sheets which have been furnished to the receiving Office in 
response to an invitation under Article 14 are referred to in this report as "originally filed" and are not annexed to 
the report since they do not contain amendments (Rules 70. 16 and 70. 17).): 

Description, pages: 

1-2 as originally filed 

Ciaims, No.: 

1 -9 as originally filed 

Drawings, sheets: 

1/3-3/3 as originally filed 

2. With regard to the language, all the elements marked above were available or furnished to this Authority in the 
language in which the international application was filed, unless otherwise indicated under this item. 

These elements were available or furnished to this Authority in the following language: , which is: 

□ the language of a translation furnished for the purposes of the international search (under Rule 23.1 (b)). 

□ the language of publication of the international application (under Rule 48.3(b)). 

□ the language of a translation furnished for the purposes of international preliminary examination (under Rule 
55.2 and/or 55.3). 

3. With regard to any nucleotide and/or amino acid sequence disclosed in the international application, the 
international preliminary examination was carried out on the basis of the sequence listing: 

□ contained in the international application in written form. 

□ filed together with the international application in computer readable form. 

□ furnished subsequently to this Authority in written form. 

□ fumished subsequently to this Authority in computer readable form. 

□ The statement that the subsequently furnished written sequence listing does not go beyond the disclosure in 
the international application as filed has been furnished. 

□ The statement that the information recorded in computer readable form is identical to the written sequence 
listing has been fumished. 

4. The amendments have resulted in the cancellation of: 

□ the description, pages: 

□ the claims, Nos.: 
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□ the drawings, sheets: 

5. □ This report has been established as if (some of) the amendments had not been made, since they have been 
considered to go beyond the disclosure as filed (Rule 70.2(c)): 

(Any replacement sheet containing such amendments must be referred to under item 1 and annexed to Uiis 
report.) 



6. Additional observations, if necessary: 



V. Reasoned statement under Article 35(2) with regard to novelty, inventive step or industrial applicability; 
citations and explanations supporting such statement 

1. Statement 

Novelty (N) Yes: Claims 1-9 

No: Claims 

Inventive step (IS) Yes: Claims 1-9 

No: Claims 

Industrial applicability (lA) Yes: Claims 1-9 

No: Claims 



2. Citations and explanations 
see separate sheet 



VII. Certain defects in the international application 

The following defects in the form or contents of the intematlonal application have been noted: 
see separate sheet 
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Reference is made to the following documents: 

D1 : NOBUG MORITA: 'Field and Laboratory Verifications of Sand Production 
Prediction Models' SPE, vol. 27341, 7 - 10 February 1994, pages 19-28, 
XP0021 46574 

D2: ZHANG ET AL.: 'Mechanical Strength of Reservoir Materials: Key 

Information for Sand Prediction' SPE, no. 49134, 27 - 30 September 1998, 
pages 423-430, XP000937927 

Re Item V 

Reasoned statement under Article 35(2) with regard to novelty, inventive step or 
industrial applicability; citations and explanations supporting such statement 

V-1 D1 , which is considered the closest prior art, discloses the subject-matter of 

claim 1 as follows: 

A method of predicting the failure of a rock formation surrounding a 
subterranean cavity (p. 19, last section I. 3-5 and last 4 lines), comprising the 
steps of 

-measuring a set of parameters relating to pressure conditions and stresses 
in the rock formation surrounding the subterranean cavity (p. 24, last section, 
lines 1-3 and part (c) and (d)); 

-using the set of parameters to determine a rock strength (p. 24, last section, 
lines 1-3 ); 

-determining a first characteristic length relating to the size of the cavity (p. 
22, left column, see text under (1)); 

-determining a second characteristic length relating to the grain size of the 
rock formation surrounding the cavity (p. 22, left column, see text under (1)); 
-correcting said rock strength (page 24, right column, I. 5-6; and 
-using a failure criterion* and the corrected rock strength to predict a 
condition under which the rock formation is expected to produce debris 
(page 19, right col., first eight lines of second paragraph). 

*) According to the applicant, the use of a failure criterion is already known, 
see the description on page 2, 1. 27-30. 
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The apparatus according to clainn 1 , therefore differs with respect to D1 in 
using the first and second characteristic lengths to determine a correction for 
the rock strength. 

The apparatus according to claim 1 is therefore new and the claim meets the 
novelty requirements of Article 33(2) PCT. 

V-2 The distinguishing features of claim 1 relates to the underlying problem of 
obtaining an indicator of relevant downhole physical dimensions for error 
correction. 

The fact that the determination of rock strength is prone to error is clear from 
both D1 (p. 25, Conclusions) and D2 (p. 424, left col. first section, last 8 
lines). 

In particular, D1 teaches that such errors should be corrected (p. 24, right 
column, I. 1-6). But even when the skilled man considers what is said on 
page 22, left column, step (2) "cavity size is not close to grain size", it would 
not be obvious to him to take these two factors and use them for error 
correction. Neither of the other cited documents in the International Search 
Report provide hints in this direction and therefore the method step of claim 
1 of using said characteristic lengths as correction involves an inventive step 
such that claim 1 meets the requirements of Article 33(3) PCT. 

V-3 Since claims 2-9 depend upon claim 1 , these claims also meet the 

requirements of Article 33(2) and (3) PCT. 

Re Item VII 

Certain defects in the international application 

VII-1 Independent claim 1 is not in the two-part form in accordance with Rule 

6.3(b) PCT, which in the present case would be appropriate, with those fea- 
tures known in combination from D1 being placed in the preamble (Rule 
6.3(b)(i) PCT) and with the remaining features being included in the 
characterising part (Rule 6.3(b)(ii) PCT). 
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If, however, the applicant is of the opinion that the two-part form would be 
inappropriate, then reasons therefor should be provided in the letter of reply. 
In addition, the applicant should ensure that it is clear from the description 
which features of the subject-matter of claim 1 are already known in 
combination from D1 (see the PCT Guidelines. Ill-2.3a). 

VII-2 Contrary to the requirements of Rule 5.1 (a)(ii) PCT, the relevant background 
art disclosed in D1 and D2 are not mentioned in the description, nor are 
these documents identified therein. 
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lol.ntifKj to prcFisure conditions and 

stresses in the rock formation 
surrounding the subterranean cavity 



determine rocic strength using the 
set of parameter 



determining a first characteristic 
length relating to the size of the 
cavity 



determining a second 
characteristic length relating to 
the grain size of the rock 
formation surrounding the cavity 



using the first and second 
characteristic lengths to 
determine a correction for the 
rock strength 



correcting said rock strength 



using a failure criterion and the 
corrected rock strength to predict 
a condition under which the rock 
formation is expected Co produce 
debris 



(57) Abstract: A method of predicting the failure of a rock forma- 
tion surrounding a subterranean cavity, including measuring a set of 
parameters relating to pressure conditions and stresses in the rock 
formation surrounding the cavity; using the set of parameters to de- 
termine a rock strength; determining a first characteristic length re- 
lating to the size of the cavity; determining a second characteristic 
length relating to the grain size of the rock formation surrounding the 
cavity; using the first and second characteristic lengths to determine 
a correction for the rock strength; correcting said rock strength; and 
using a failure criterion and the corrected rock strength to predict 
a condition under which the rock formation is expected to produce 
debris. The results of the prediction can be used to monitor wellbore 
stability while drilling or optimize the production parameters for a 
hydrocarbon reservoir. 
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Cavity Stability Prediction Method for Wellbores 



This invention relates to a method of estimating or predicting 
the stability of cavities in a subterranean formation. It 
5 further pertains to using such estimates to control and set 
operation parameters for drilling and producing hydrocarbon 
wells . 

10 BACKGROUND OF THE INVENTION 

For the production of hydrocarbon wellbores are drilled into 
subterranean formations. Subsurface formations encountered in 
oil and gas drilling are compacted under in situ stresses due to 

15 overburden weight, tectonic effects, confinement and pore 

pressure. When the wellbore is drilled in a formation, the rock 
near the wellbore is subjected to increased shear stresses due 
to a reduction in confinement at the wellbore face after removal 
of the rock from the hole. Compressive failure of the rock near 

20 the wellbore will occur if the rock does not have sufficient 
strength to support the increased shear stresses imposed upon 
it. 

Formation stability problems are not only encountered during the 
25 drilling of the wellbore. For the production of hydrocarbons, 
the hydrocarbon bearing formation is usually perforated or 
fractured to enable and stimulate the fluid flow into the 
wellbore. When producing from unconsolidated or weakly- 
consolidated reservoirs, the formation tends to produce 
30 particulates, (e.g. sand) along with the hydrocarbons. 



Formation sand is produced when the combined effects of fluid 
drag and near-wellbore stresses cause disaggregation near the 
perforation or fracture. Individual grains of sand are detached 
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from the matrix forming the formation. At relatively low flow 

rates, fluid drag does not affect the stability, but as flow 

rate increases, drag forces become sufficiently high to remove 
sand particles from the matrix. 



Flowrate from a formation is normally controlled by the 
perforation drawdown pressure (DP) which is the difference 
between the pore pressure (p^) in the formation and the 
bottomhole pressure (Pq) and can hence be expressed as DP= Po - 



The critical drawdown pressure (CDP) is the value of DP at which 
the rock matrix surrounding the perforation begins to de- 
stabilize. Its value is determined by the maximum calculated 
15 rock strength. 

To model the maximum rock strength classical elastic and elasto- 
plastic theories, failure criteria and fracture mechanics have 
been applied. Models use empirically or semi-empirically derived 
20 rock strength values to predict formation behavior by using 
classical theories and stress, pore pressure and empirically 
derived strength data from various wells . 

There are several methods for predicting when for example sand 
25 production will occur in a particular well. Such methods are 
disclosed and discussed in the US Patent No 5,497,658 and 
references contained therein. Known rock failure criteria as 
discussed in this and other published docoiment are referred to 
as Mohr-Coulomb, critical state, Drucker-Pager model or as 
30 extended Von Mises criterion 

To apply the failure criteria it is necessary to measure rock 
properties and the formation fluid properties from core samples, 
wellbore logs, and the like. 



5 



10 pw. 



35 
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It is therefore an object of the present invention to provide a 
novel method of estimating the strength of cavities in the 
subterranean formation, particulary the initiation of sand 
production in subterranean (sandstone) formations. 



According to one aspect of the present invention, there is 
10 provided a method of predicting the failure of a rock formation 
surrounding a subterranean cavity, including the steps of 
measuring a set of parameters relating to pressure conditions 
and stresses in the rock formation surrounding the cavity; using 
the set of parameters to determine a rock strength; determining 
15 a first characteristic length relating to the size of the 

cavity; determining a second characteristic length relating to 
the grain size of the rock formation surrounding the cavity; 
using the first and second characteristic lengths to determine a 
correction for the rock strength; correcting said rock strength; 
20 and using a failure criterion and the corrected rock strength to 
predict a condition under which the rock formation is expected 
to produce debris . 

A cavity can be a wellbore without lining (open hole) or 
25 perforation tunnels or other spaces created in a subterranean 
formation by using chemical or physical forces such as 
explosives and drilling equipment. 

The set of parameters used to characterize the formation 
30 surrounding the cavity may include measurement as performed by 
logging devices, such as sonic, gamma-ray logging devices or NMR 
based logging devices. Important parameters are for example 
density or porosity, clay content, or p- and s-wave slowness. 



5 



SUMMARY OF THE INVENTION 
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The characteristic length relates to the dimensions of a cavity 
or grain and is preferably the diameter or radius or the closest 
approximation of the diameter or radius, given the irregular 
dimensions of those subterranean objects. 

The results of the prediction can be used to monitor wellbore 
stability while drilling or optimize the production parameters 
for a hydrocarbon reservoir. 

The normalization of the cavity dimension or length with the 
grain size yields a correction factor that can be used to derive 
an apparent rock strength- In this way, the scale and plasticity 
effects are lumped into an apparent strength calculation. This 
apparent rock strength can be used with estimates of in- situ 
stresses and pore pressure in a 3-D poroelastic model and 
failure criterion as Mohr-Coulomb for the calculation of the 
critical parameters related to the stability of the cavity, such 
as draw-down pressure and the onset of sand production. 

Combined with the appropriate measuring-while-drilling (MWD) or 
logging-while-drilling (LWD) technology, it can be converted 
into a prediction tool to estimate the rock stability during 
drilling operation in real time. As such it could contribute 
significantly to the prevention of s tuck-pipe problems, 
currently the cause of significant losses in the oilfield 
industiry . 

These and other features of the invention, preferred embodiments 
and variants thereof, possible applications and advantages will 
30 become appreciated and understood by those skilled in the art 
from the following detailed description and drawings. 



35 
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BRIEF DESCRIPTION OF THE DRAWINGS 



5 



FIG. 1 



is a schematic drawing of a wellbore and a perforation 
tunnel illustrating the directions of stresses; 



FIG. 2 



shows the critical draw-down pressure curve for a 
simulated reservoir; and 



10 FIG. 3 



charts steps of the present invention. 



MODE(S) FOR CARRYING OUT THE INVENTION 



The underlying idea is to use log-data (mainly sonic data) for 
15 the derivation of rock elastic constants and formation strength 
parameters. These parameters can be used with estimates of in- 
situ stresses and pore pressure in a 3~D poro-elastic model and 
Mohr-Coulomb failure criterion for the calculation of the 
critical draw-down pressure. 



The method described below assumes clean sandstone as formation 
material . 

The bulk porosity can be derived from the bulk density pb of a 
25 fluid saturated porous rock, which is given by 



where Ps is the density of the solid grains and Pf .is the fluid 
30 density. Solving for the bulk porosity results in 



20 



[1] 



Pb = 9Pf + (1 - <P)Ps 



(p = - Pb 
Ps - pf 
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Approximate default values can be assumed for both densities, 
e.g., Ps = 2.75 g/cm^ and Pf = 1.1 g/cw?. 

The elastic parameters are computed from log compressional and 
5 shear wave velocities. Methods and apparatus to perform the 

required measurements are known as such in the art. For example, 
the United States Patents 4,862,991, 4,881,208 and 4,951,267 
refer to logging tools for measuring shear and compressional 
wave slowness. The Schlumberger DSI™ tool for conventional 
10 logging or the ISONIC™ tool for logging-while-drilling are 

capable of measuring the required data. Reference to those tools 
are found for example in the Schlumberger Oilfield Review, 
Spring 1998, 40-66. 

15 The elastic parameters of the formation as used by the present 
invention can be determined using the compressional and shear 
wave velocities log data. The Poisson ratio v, the shear modulus 
G, the Young's modulus E and the bulk modulus K are calculated 
from the p and s wave slownesses (i.e. the reciprocal of the 

20 velocity) , Dtc and Dts, according to equations: 



[3] 



0,5(Dt3 / Dt^)^ - 1 
(Dt3 /Dt^f - 1 



[4] 



Pb 




25 



[5] 



2G(1 + V) 



[6] 



3(1 - 2V) 
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The rock strength parameters can be calculated in terms of the 
uniaxial (or unconfined) compresssive strength UCS from the 
empirical correlations known as Coates and Denoo equation: 

5 [7] UCS = (114 + 97Vsh) K(in mio. psi) E(in mio . psi) 

where the clay content Vsh can be determined using for example 
gamma ray logs or information from core. 

10 The pore pressure, Pq, is given by the reservoir pressure . 

Methods and apparatus to measure the reservoir pressure (and the 
wellbore pressure Pw ) are known and reference is made to the 
United States Patent 5,789,669 for details of such measurements. 
The reservoir pressure is likely to vary with time according to 

15 the predicted performance of the reservoir. 

The vertical in-situ stress (illustrated by FIG. 1) is 
estimated from the overburden weight. The magnitude of the 
minimiim horizontal stress can be obtain either from 
20 consolidation theory according to 




where p is the Biot coefficient, or from frictional equilibrixim. 
25 If possible, a stress measurement or extended leak-off test 
should be used to verify which assumption gives better 
estimates . 

Finally, in a tectonic environment the horizontal stresses are 
30 unequal 



a„ = K Oh 
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The ratio between horizontal stresses can be estimated from 
borehole breakouts or by the simulation of field tectonic 
movement using finite elements. In general as much information 
as possible should be used in constraining the values of the 
horizontal stresses • 

In the following the methodology for calculating the optimiim 
draw-down pressure DP based on 3-D elastic solution. The basic 
equations are known. The known 3-D elastic solution is augmented 
with extra terms for taking into account for the gradient of 
pore or reservoir pressure during production. 

As illustrated by FIG.l, the method can be applied to estimate 
the stability of sections of the wellbore or to estimating the 
stability of other cavities such as perforation tunnels . 

Transforming the parameters from a vertical into a wellbore 
coordinate system, the stresses at a point on the borehole wall 
(r = R) and at an angle 0 from the axis x are given by 

tlO] = 

C^e = ( <yxx + O^yy - Pw) " 2( - Gyy) COS 29 - 

« 1 - 2v 

- 4a^ sin 29 - (Po - p„)P -^—^ 
f^z = <^zz - 2v( a,^ - Oyy) COS 29 - 

[12] 1 - 2V 

- sin 29 - (Pq - p„)P 
[13] Obz =- 2a^2 sin 9 - 2a„, cos 9 
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[14] G^^ = 0 

where the original input in-situ stresses, Qh, g^, cJv have first 
been transformed into the Cartesian components of a wellbore 
coordinate system and then, using eqs [10]- [14], into 
cylindrical wellbore coordinates. The parameter p„ denotes the 
pressure in the wellbore. For a wealc reservoir sandstone a 
reasonable value for the Biot coefficient is P = 1. 

The principal stresses can be found from the eigenvalues of the 
stress tensor 











^rz 


[15] 


[a] = 






Oez 










€^z. 



using the Matlab™ function princ = eigs(s), and can be put in 
order, G^, O2 and ai,the maximum compressive stress. 

The Mohr-Coloumb failure criterion can be expressed in the 
following form 

[16] f = UCS - G\ 

The effective stress G'l at the borehole wall is given by 
[17] G\ = - Pp^ . 

It was found that the failure criterion, eq. [16] , and any other 
failure criterion using the uniaxial compressive strength UCS 
can be improved by talcing into account the scaling effect, i.e. 
the characteristic dimension of the perforations through which 
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hydrocarbons are produced. Experimental data showed that by 
introducing a scaling factor including the grain size of the 
formation, the estimates of the critical production parameters 
can be improved and applied to a broader range of rock types. 

Applying the scaling factor to the uniaxial compressive strength 
UCS yields the correction 



where UCS is defined by eq. [7] and Dperf is the diameter of the 
perforation and Dgrain is the diameter of the grains of the rock 
foirmation. The fitting parameters a and n are determined as 
16.1064 and 0.3374, respectively, by may vary to some extend 
depending on the fitted data and fitting algorithm. 

In the absence of a measured grain size, Dgrain can be estimated 
using prior knowledge of the rock or, at worst, simply 
approximated by a constant default value. Experimental data 
suggest 0.2 mm for such a default value. 

The corrected UCSappar. can be used in the failure criterion [16] 
and standard mathematical optimization procedures to produce a 
better estimate of the maximal rock strength and, hence, a 
better estimate of the maximum draw-down pressure. 

FIG 2 illustrates a simulated example using input values taken 
from known parameters of a drilled well in the North Sea. 

The input parameters are 




-n 



[18] 



UCS 



'appar. 



= 2 UCS a 



Insitu stresses: 
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Vertical stress Gv = 24.82 MPa; 
Min. horizontal stress = 15.63 MPa; 
Max. horizontal stress Gh = 17.19 MPa; 
Formation pressure Pq = 11.03 MPa. 

Rock Parameters: 
Poisson ratio v= 0.25; 

Uniaxial compressive strength UCS = 4.07 MPa; 
Grain size 

J-'grain — U . ^ mm 

Well data: 

Well diameter Dweii = 0.20 m 
Inclination I = 90 degrees 
Azimuth a= 0 degrees 

Perforation data 

Perforation diameter Dperf = 0.01 m 
Phasing (()= 55 degrees 

The horizontal stresses are assumed to be equal and they are 
calculated from the consolidation eq. [9] . The formation 
strength is calculated in terms of the corrected UCSappar. from 
available log data and the correlation function [7] • 

FIG, 2 shows the optimum wellbore pressure for sand-free 
production calculated using the above approach at the beginning 
of (0% depletion) and during production. During depletion it is 
assumed that the total vertical in-situ stress remains 
unchanged, therefore, the vertical effective stress increases by 
the same amount the pore pressure decreases. The variation of 
the effective horizontal stresses is taken empirically to be 50% 
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of the variation in the vertical effective stress. Though safe 
production is possible within the area limited by calculated 
curve for the onset of sand production (marked by circles) , 
maximum hydrocarbon is achieved by setting the well parameters, 
5- i.e. most notably the wellbore pressure as close to the curve as 
possible . 

Using the same input data and stability model (i.e. UCS) without 
the correction proposed by the present invention, the 
10 optimization predicts that the wellbore can not be produced 
without sand- 
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CLAIMS 



10 
15 
20 
25 

2. 

30 

3 . 



A method of predicting the failure of a rock formation 
surrounding a subterranean cavity, comprising the steps of 

- measuring a set of parameters relating to pressure 
conditions and stresses in the rock formation surrounding 
the subterranean cavity; 

- using the set of parameters to determine a rock strength; 

- determining a first characteristic length relating to the 
size of the cavity; 

- determining a second characteristic length relating to 
the grain size of the rock formation surrounding the 
cavity; 

- using the first and second characteristic lengths to 
determine a correction for the rock strength; 

- correcting said rock strength; and 

- using a failure criterion and the corrected rock strength 
to predict a condition under which the rock formation is 
expected to produce debris- 

The method according to claim 1 wherein the set parameters 
includes sonic wave slowness. 

The method according to claim 1 wherein the set parameters 
includes the formation density. 



4. 

35 



The method according to claim 1 wherein the set parameters 
includes the wellbore and formation pressure. 
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5. The method according to claim 1 wherein the failure 

criterion is a shear failure criterion (Mohr-Coulomb) . 

5 6. The method according to claim 1 wherein the failure 

criterion includes a term corresponding to an uniaxial 
compressive strength (UCS) . 

7 The method according to claim 1 wherein the correction 
10 includes forming the c[UOtient of the first and the second 

characteristic length , 

8. The method according to claim 1 further including the step 
of determining a wellbore production pressure using the 

15 failure criterion. 

9. The method of claim 1 wherein the set of parameters 
relating to pressure conditions and stresses in the rock 
formation surrounding the cavity are at least partly 

20 measured while drilling. 
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Measuring a set of parameters 
relating to pressure conditions and 
stresses in the rock formation 
surrounding the subterranean cavity 



determine rock strength using the 
set of parameter 



determining a first characteristic 
length relating to the size of the 
cavity 



determining a second 
characteristic length relating to 
the grain size of the rock 
formation surrounding the cavity 



using the first and second 
characteristic lengths to 
determine a correction for the 
rock strength 







correcting said rock strength 







using a failure criterion and the 



corrected rock strength to predict 
a condition under which the rock 
formation is expected to produce 
debris 



FIG. 3 
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